UBIVERBHT OF TTLTNOIS 
DIGITAL COMPUTER 



LIBRARY ROUTB E 5 5 - 85 



TITLE 

TYPE 

HIMBSR OF WORDS 

TEMPORARY STORAGE 



ACCURACY 
DURATION 



PRH5ET PARAMETER 



Minimization of a Function of n Variables (DOI or SADOl) 
Closed 

89 

n + 1 words at S3 

n + 1 words at Sk 

n + 1 words at S5 

Depends upon the condition of the function 

A "minor" cy«le takes T = («5 + l^n 4- t) us where n is the 

number of variables and t = time to compute f(x^, «.., x^)* 

A "major" cycle takes (1.8 + h.Jn + 2nt + HT) where N is 

the numb*r of "minor cycles" in a major cycle. 

S3 - SH 

where a, a+1, .-., a+n are the addresses 



3 - 00 F 00 aF 

k - 00 F 00 bF 

5 * 00 F 00 cF 

6 - 00 F 00 OJ 

7 - 00 F 00 eJ 

8 - 00 F 00 sF 



C 2* 



x on 

n 



of f\X- f •••> x ) f x^i 
entry and exit $ 
b, b+1, ..., b+x are the addresses of 






5 * 
*n 



c, e+1, w»-, c+n are working spaces j 

a is the factor by which the mesh size 

is decreased; 

e is an end constant such that only mesh 

siaes larger than € will be usedj 

s is the address of a closed subroutine 



which takes x 1 from c+1 (1S5), Xg from 
c+2 (2S5), . .., x from c+n (nS5) and 
places f(x L , Xg, ..., x q ) in R.^ 

9 - 00 F 00 (a+n)F n is the number of -variables 
K - 00 F 00 (b+n)F 
S - 00 F 00 (c+n)F 

Iff - 00 F 00 rF or The preset parameter in N allows the 
00 F ^0 F programmer some choice in scalej if 



-2- 



It is 00 F 00 rF 

5 VILLI "be computed asi 

X i 

2- r f 

8 ~ > i i i i ******** 

1 S (2~ r f 

I X I 



whereas if it is 00 F 40 F we set r = in the ahove 
expression. In the ahove* 

f = f(.««x. ■» B«*«) ** f(«.«x. + 5«.«) 

In either case if the denominator exceeds scale it is 

-59 
replaced hy 1 - 2 . 

ERFBX x^ x , ... x^ the initial values of x*, X-j ..., x 

are to "be loaded in 1S3, 2S3, ••-, nS3. 8 , the initial 

value of 5 is to he in IL , and the rain program should 

contain: 

50 rF 

50 pF 



Jh-1 



26 qF 



where q is the address of this routine and r is the address 
to the left hand side of which control will he transferred 
"before every decrease in mesh size. At this address the 
programmer may place a routine to assess or print inter- 
mediate results. This routine should return control to 
the left hand side of (qVjp-)* This can he done automatically 
since when control is transferred to r(L.B. ) the appropriate 
return address will he in the right hand address position 
of R 2 and can he utilized hy S5, 42 into a 26 at the end" 
of the programmer's interlude routine. The hast c urr en t 
values of ffx^ - .., x^), x^, .«., x^ are always in S3, 
1S3, • ««, nS3« If it is not desired to leave this code, 
set r = q + 71. The routine Is finally left with the 
hest values of ffe^ Xg, ..., x^), x^, x^, ..., x in 
S3, 1S3, 2S3, ..., nS3 respectively. 
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BESCRIFriGir OF METHOD 1 - f are computed 

x i 

2-5 are computed 

x i 



HOPE 



3 - f(x- +5 , ..., x+ 6 )is computed and tested 
1 a 

against previous ffo* •••* x q ). If the function 
hym decreased the old 'values of function and argument 
are replaced "by the new* The process is repeated 
with the same 5 . This step comprises a "minor" 
cycle . 

k - When the above process ceases to improve the function, 

we return to step 1 and compute new 6 . Steps 

1 
1-4 comprise a "major" cycle. 

5 - When the first "minor" cycle of a given "major" 

'cycle fails to Improve the function, indicating 
that no further Improvement is to "be expected 
using the current mesh siae, control is transferred 
to the interlude as explained above. 

6 - Upon returning to the code, 8, the mesh size will be 

replaced by aS -where a is entered in S6 during read- 
in. The steps 1-5 are then repeated. 

7 - Steps 1-6 will he repeated (54-1) times where U is 

the largest integer for which a 5 Q > €• € is to 

he entered in S7 during read-in. 
A function, poorly conditioned in having a very small gradient 
with respect to some argument, can deceive this routine 
in that this coordinate of the minimum will he very poorly 
found. If the programmer suspects that such a condition 
exists, this routine can he used for roughing and Routine 
H 6 - 86, a brute force approach to the minimum by varying 
one argument at a time, may be used for finishing. All 
parameters (except SK and SN), entry and contents of 
storage location of this code and Routine H 6 are identical, 
so that they may be used at the same time in the machine 
for the above purpose* 
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Routines H 3, H k, H 5, and H 6 constitute a set sufficient 
(perhaps not necessary) to minimize any function of n 
variables. Dae to economy of storage, Routine H 3 should 
he used for n = 1 or 2, Routine H k for n = 3 or k, and 
Routine H 5 for n > 5« 



lgr 
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LOCKTXOH 



3 
4 

5 
6 

7 

8 

9 

10 
11 
12 

13 

14 

15 
16 



ORDER 



00 K(H5) 
1*0 84 

S5 IF 

46 60L 
L4 L 
42 73L 
L5 3L 
50 S8 

46 10L 
L5 77L 
40 5L 
L5 F 
40 T 
L5 5L 

l4 j6l 

40 5L 

L5 80L 
LO 5L 
36 5L 
50 131 
50 9L 
26 F 
40 S3 
L5 9L 

47 10L 
26 4L 
00 F 
L5 60L 
46 6% 
L5 77L 
46 19L 
46 23L 
40 27L 
L5 22L 
46 20L 



KJEB 



BUSK 1 



1 5 



From 12, 73 

From 8 
:$r 4, 7 



By 3, H 



From 10 



Store 5- 

Moff» lit link 

Plant interlude address 

Plant 11 nX 

Set 1C3L for initial entry 



Set initial addresses in transfer orders 



x£, ,,«, x£ to 1S5, «*•, »S5 



j^ranse addresses in transfer orders 



Teat far i = n 



f(3c_, .•., x n ) to S3 



Set 10L for subsequent entry 



Set snitch for failure on first try 



Set i « 1 in 19 - 27 



LOCATION 
IT 

18 

19 
20 
21 
22 

23 
2k 

25 
26 

27 
28 

29 
30 
31 
32 
33 



GEEER 

1*6 2l*L 
L5 T9L 
1*2 22L 
1*0 26L 
L5 F 
lU s4 
1*0 F 
50 20L 
26 S8 
10 SB 
50 1S5 
1*0 F 
L5 F 
LO Si* 
1*0 F 
50 2l*L 
26 S8 
10 SH 
LO F 
1*0 F 
L5 F 
1*0 F 
L5 2TL 
Ll* T6L 
1*0 27L 
h6 19L 
h6 23L 
L5 20L 
Ll* T6L 
k6 20L 
k6 2l*L 
L5 26L 
L5 76L 
1*0 26L 



noras 



BIGE 2 



H 5 



By 1*, 29 
From 55 
By 16, 31 



By 18, 3^ 
By 15, 30 

By 17, 32 



By 18, 33 

By 18, 33 

By 15, 29 

By 15 , 29 



Ccnprccte T 



Achnance i to i+1 in 19-27 



LOCATION 


ORDER 




HOTES BW5E 3 


3* 


1*2 22L 








L5 80L 




Test far I = n 


35 


LO 27L 
36 19L 






36 


L5 79L 








46 38L 




Set i = I in"38 


37 


19 38F 
40 S5 1 






38 


L7 F 
L4 S5 


By 36, 40 
From 42 




39 


40 S5 
L5 38L 






40 


L4 76L 




Ccnrptrte £ |f J 
i x i 




46 38L 




in 


L5 81L 
LO 38L 






Ha 


36 38L 
26 8% 






*3 


46 44L 
42 4% 


From 87 


Set i* = 1 in 44-45 


J 44 


L5 F 


3y 43, 47 






66 S5 




Caarptrte 6 

x l 


^5 


7J S4 






40 F 


By 43, 48 




46 


L5 44l 






; 


L4 76L 




Advance i to i+1 tn 44-45 


! 


46 44l 
L5 4% 






J 48 


L4 T6L 
k£ 45L 






! 49 


L5 82L 






i 


LO 4% 




Test for i = n 


50 


36 44L 






1 


L5 77L 

4 


From 70 





H 5 











wcssim 1 


ORDER 


ions nss 4 


51 


46 5to 
L5 T9L 




Set 1 - 1 ia 5**«"55 


52 

! 


42 5^ 
L5 22L 




1 

1 


53 


46 55L 

50 7IL 




j 


5* 


L5 F 


By 51, 56 

From 59 


• 




l4 f 


By 52, 56 


Set x 4 + b la 1S5 
1 *i 


55 


40 F 


By 53, 58 




L5 54L 






56 


L4 76L 








40 54l 




Advance i to 1+1 ia 54-55 


5? 


L5 55L 
L4 76L 


' 




58 


46 55L 








L5 85L 




j Teat for i = a 


59 


LO 55L 
36 54L 




! 


60 


50 F 


By 1 


Interlude address 




50 6ql 




f(.,.x. + 5 .*.) 
1 x i 


61 


26 S8 




3 

I 


J 40 S5 




Form f *• f (before) 


| 62 


LO S3 
50 53L 






63 


36 F 


By 13, 69 


To iaterlade 




L5 78L 




* 


64 


42 65L 
46 65L 




Set I * 1 la 65 


65 


L5 F 


By 64, 67 
From 68 






40 F j 


By 64, 67 


Replace old valaes oy aev 


66 


L5 6% 








l4 76L 




Adraace 1 to 1+1 la 65 



H 5 



mmiGR 


ORDER 


61 


ko 65L 




L5 &L 


68 


LO 65L 




36 6^L 


69 


L5 7QL 




k6 63L 


70 


50 !4l 




22 5QL 


71 


50 T^L 




7J S4 


72 


ko Sk 




LO 75L 


73 


36 kL 




22 F 


7^ 


00 F 




00 S6 


75 


00 F 




00 S7 


76 


00 IF 




00 IF 


77 


L5 1S3 




ko 1S5 


78 


00 S5 




00 S3 


79 


LO 1S4 


? 


ko lSk 


80 


L5 S9 




4o SS 


81 


L7 SK 




L4 S5 


82 


7J S4 




40 SK 


83 


40 SS 




L5 5^L 



IDTES 



KAfig 5 



I 5 



Test for i ■ a 



Open fttilare Jralteh and repeat 



By 2. 



Minor cycle 



Foam 8 



141 -^ 



fkn& test against e 



Link 



= a 



= € 

Address advance 



Starting constants 



End. constants 



i LOCATION 


ORDER 




8k 


L5 SS 
kO S9 




85 


L5 S5 
56 87L 


From ^2 


86 


L5 88L 
kO S5 




87 


L5 79L 
26 43L 




88 

1 

1 

1 

i 
J 

! 
} 

i 

> 

i 


7L 4095F 
LL 4095F 

> 


L 



ions 



PkCEE 6 



H 5 



Test S If 



I X i 



If > 1, replace fcy 1-2 
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J 



